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Abstract. To determine nonspherical angular momentum amplitudes in hadrons at long ranges (low Q^), 
data were taken for the p(e, e'p)7r'' reaction in the A region at Q^ — 0.060 (GeV/c)^ utilizing the magnetic 
spectrometers of the Al Collaboration at MAMI. The results for the dominant transition magnetic dipole 

3/2 

amplitude and the quadrupole to dipole ratios at W — 1232 MeV are: M^^ — (40.33 ± 0.63stat+syst ± 

0.6Uodel)(10"V"^^+), Re(Si{.VA^ff ) = (-2-28 ± 0.29,tat+.y.t ± 0.20,nodei)%, and Re{Sli^/Mfi^) = 
(—4.81 ± 0.27stat+syst ± 0.26modoi)%. These disagree with predictions of constituent quark models but are 
in reasonable agreement with lattice calculations with non-linear (chiral) pion mass extrapolations, with 
chiral effective field theory, and with dynamical models with pion cloud effects. These results confirm the 
dominance, and general Q^ variation, of the pionic contribution at large distances. 

PACS. 1 3.60.Le, 13.40. Gp, 14.20. Gk 

Experimental confirmation of the presence of nonspher- duction experiments at JLab [6l[7| for Q^ from 0.4 to 4.0 

ical hadron amplitudes (i.e. d states in quark models or (GeV/c)^, at Bates [8ll9 l ll0 p if] at Q"^ = 0.127 (GeV/c)^ 

p wave TT-N states) is fundamental and has been the sub- and Mainz [T21ll3l[Tl] at Q^ = 0.127, 0.20 (GeV/c)^ have 

ject of intense experimental and theoretical interest (for been published. The present result at Q^ = 0.060 (GeV/c)^ 

reviews see [UEllS]). This effort has focused on the mea- is the lowest Q^ value probed to date in modern, preci- 

surement of the electric and Coulomb quadrupole ampli- sion electroproduction. It adds a very important point to 

tudes (E2, C2) in the predominantly Ml (magnetic dipole- determine the physical basis of long range nucleon and 

quark spin flip) 7*iV — > Z\ transition. Measurements of the A nonspherical amplitudes and is a test of the Q^ region 

E2 amplitude from photopion reaction experiments with where pionic effects are predicted to be dominant and ap- 

polarized photons have been reported [HIS] . Electropro- preciably changing. 

Since the proton has spin 1/2, no quadrupole moment 
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tion multipolcs are M^\^ , E^j^ , and S^j^ . The relative 

0/9 0/0 

quadrupole to dipole ratios are EMR = Rc(i?]^^ /-^'^r+ ) 

and CMR = Rc(S'i|V*^iT)- In the quark model, the non- 
spherical amplitudes in the nucleon and A are caused by 
the non-central, tensor interaction between quarks [15] . 
However, the magnitudes of this effect for the predicted 
E2 and C2 amplitudes [1^] are at least an order of mag- 
nitude too small to explain the experimental results (see 
Fig.[3]below) and even the dominant Ml matrix element is 
~ 30% low inUm ■ A likely cause of these dynamical short- 
comings is that the quark model does not respect chiral 
symmetry, whose spontaneous breaking leads to strong 
emission of virtual pions (Nambu-Goldstone Bosons) [2. 
These couple to nucleons as 1? • p where a is the nucleon 
spin, and p is the pion momentum. The coupling is strong 
in the p wave and mixes in non-zero angular momentum 
components. Based on this, it is physically reasonable to 
expect that the pionic contributions increase the Ml and 
dominate the E2 and C2 transition matrix elements in 
the low Q^ (large distance) domain. This was first indi- 
cated by adding pionic effects to quark models [T7], sub- 
sequently shown in pion cloud model calculations [I81IT9] . 
and recently demonstrated in chiral effective field theory 
calculations POU2T] . 

The five-fold differential cross section for the p(e, e'p)7r° 
reaction is written as five two-fold differential cross sec- 
tions with an explicit (jf dependence as |22| 



rfV 



dQfdEfdQ 



r{(JT + f^f^L + Vlt<7LT cos (j)* 

eaTT cos 20* -f hpeVLT'<7LT' sin (j)*) (1) 



where e is the transverse polarization of the virtual pho- 
ton, vlt = v^2e(l 4- e), vlt' = A/2e(l — e), F is the vir- 
tual photon flux, 0* is the pion center of mass azimuthal 
angle with respect to the electron scattering plane, h is the 
electron helicity, and Pe is the magnitude of the electron 
longitudinal polarization. The virtual photon differential 
cross sections {ar, cl, ctlt, o'tt, o'lt') are all functions of 
the center of mass energy W , the four momentum transfer 
squared Q^, and the pion center of mass polar angle 0* 
(measured from the momentum transfer direction). They 
are bilinear combinations of the multipolcs }22j . 

The p{e, e'p)^'^ measurements were performed using 
the Al spectrometers at the Mainz Microtron[23]. Elec- 
trons were detected in Spectrometer A and protons in 
Spectrometer B. Timing and missing mass cuts were suf- 
ficient to eliminate the tt" background. Spectrometer B 
has the ability to measure at up to 10 degrees out-of-plane 
in the lab. Due to the Lorentz boost this is significantly 
larger in the center of mass frame. The Mainz Microtron 
delivers a longitudinally polarized, continuous, 855 McV 
beam. Beam polarization was measured periodically with 
a M0ller polarimeter to be « 75%. The beam of up to 25 
^A was scattered from a liquid hydrogen cryogenic target. 
Sequential measurements were made at M^ = 1221 MeV, 
g2 = 0.060 GeVVc^, and 6^^ = 143, 156, 180°. For non- 
parallel measurements (6* ^ 180), the proton arm was 



moved through three (j)* settings while keeping the 0* 
value constant. The spectrometers were aligned in the lab 
with a precision of 0.6 mm and 0.1 mrad with a central 
momentum resolution of 0.01% and angular resolution at 
the target of 3 mrad [55]. The beam energy has an ab- 
solute uncertainty of ±160 keV and a spread of 30 keV 
(FWHM) [23]. The effects of these uncertainties and the 
various kinematic cuts were studied to estimate an overall 
systematic error for the cross sections of 3 to 4%. This was 
tested with elastic electron-proton scattering and the data 
agree with a fit to the world data [H] at the 3% level. 

With measurements at three 0*^ values at a fixed 0*^ 
and using a polarized electron beam, the cross sections 
(JQ ~ <Jt + ecL, (Ttt, O'lt, and aLT' can be extracted from 
the (p* and beam helicity dependence of the cross section. 
Care was taken to ensure good kinematic overlap between 
the different angular settings. Typically the phase space 
overlaps were AW « 40 MeV, AQ'^ « 0.04 GeVVc^, 
A9lg w 10°, and Acj)'^^ « 40°. Since the cross sections 
vary across the spectrometer acceptance, the shape of the 
cross section given by several models was used to refer all 
of the points to the center of the acceptance. This is a 
small correction (typically 3%) and depends only on the 
relative cross sections across the spectrometer acceptance. 
Several models were used for the collapse and each gives 
results consistent at the 0.5% level f25j . 

The measured partial cross sections are plotted in Figs.[T] 
and [2] Figure [T] shows ao and aLT with the chiral effective 
field theory (EFT) predictions [21] which have uncertain- 
ties which reflect an estimate of the neglected higher order 
terms in the chiral expansion. The other models are the 
phenomenological model MAID 2003 [28], the pion cloud 
dynamical models of Sato and Lee [18] and of DMT [19] 
(Dubna, Mainz, and Taipei), and the SAID multipole analvsis[29]. 
There is a significant spread in these model calculations 
due to differences in the resonant and background ampli- 
tudes. However it is impressive that the four model curves 
almost fall on top of each other when the three resonant 
7*p -^ A amplitudes {M^!^^ , E-^'^ , S^^ ) are varied to fit 
the data as shown in the lower panel of Fig. [1] In ad- 
dition, this panel shows the "spherical" calculated curves 
when the resonant quadrupole amplitudes (£^i{" in (Tq and 

Si_^ in O'lt) are set equal to zero. The difference between 
the spherical and full curves shows the sensitivity of these 
cross sections to the quadrupole amplitudes and demon- 
strates the basis of the present measurement. The small 
spread in the spherical curves indicates their sensitivity to 
the model dependence of the background amplitudes. 

Figure[2]shows the measured cross sections for axT and 
O'LT' with the same model curves as in Fig. [1] It is seen 
that within the relatively large estimated uncertainties the 
chiral effective field theory calculations [21] are consistent 
with experiment. For cttt the four model calculations [28) 
[TSlfTnil^ are in good agreement with each other and with 
the data. However for olt' this is not the case and only 
the Sato-Lee model agrees with experiment. 

It is instructive to examine why aLT is sensitive to 
the Coulomb quadrupole amplitude, Im[5'i+], and olt' is 
primarily sensitive to the background. The time reversal 
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even observable u^t contains the interference amplitude 
Re[5'j'^-Mi+] which is primarily sensitive to Im[S'i+]Im[Mi+ 
where the latter is the dominant multipole amplitude. By 
contrast the time reversal odd observable cflt' contains 
Im[S'J'_|_Mi+] . This is primarily sensitive to Re[S'i+]Im[Mi+] 
and therefore does not measure the Coulomb quadrupole 
amplitude but rather is sensitive to a background term 
times the dominant magnetic dipole term. The details will 
be presented in a future publication [27]. 

On the other hand, the other cross sections show a 
minimum of model dependence and can be used to ac- 
curately extract the three resonant amplitudes as will be 
discussed below. The model curves after fitting for utt 
and ult' are almost identical to those before and so have 
been suppressed, ult' is insensitive to the resonant pa- 
rameters, as mentioned above, and utt while sensitive to 
i?i+, is dominated by the Mi+ term. 

As has been discussed above, we have obtained the val- 
ues of the three resonance amplitudes {M-^i, E-^i^ , Si_^ ) 
using fits with four reaction models p8[|T8lll9il2"^ . Corre- 
lations between the fitting parameters were taken into ac- 
count in the errors estimated by the fitting routine [251 
I26j . The fits were performed using the spectrometer cross 
sections and were the same (within the errors) whether 
or not the ctlt' data were included. In addition, the fits 
used the entire / = 3/2 amplitude so that the unitarity 
of each model was preserved. At resonance, these / = 3/2 
amplitudes are purely imaginary due to the Fermi- Watson 
theorem [22]. The fits were performed at the same value 
of W at which the data were taken. The models were 
then used to extrapolate the value of the multipoles at 
W = 1232 MeV. Following our previous work[2][TT] we 
took the final results to be the average of these model de- 
terminations and estimated the model dependent error in 
the resonance amplitudes by taking the RMS deviation of 
the values [5S]. We believe that this is reasonable since 
the chosen models represent state of the art calculations 
and also a variety of different approaches. The results of 
the fits for the resonant multipoles along with the EFT 
predictions are presented in Table [T] along with the origi- 
nal values for several models. We also present the average 
fitted values for the four reaction models considered here. 
The differences between these values represent the model 
dependence due to the different background multipoles. 
The effect of background amplitudes on the resonant am- 
plitudes was studied and determined to have an effect ap- 
proximately the same size as the model to model RMS 
deviation. This study is detailed in Refs. [25] and [27]. For 
the result of this experiment we take the average values 
of the fitted multipoles using each model along with both 
the experimental and model error. It can be seen that 
the model and experimental errors are approximately the 
same magnitude. There is generally good agreement be- 
tween the EFT predictions and our experimental results. 
This also indicates the importance of the pion contribution 
to these amplitudes. It can also be seen from Table [T] that 
the dispersion between the original model calculations for 
the quadrupole amplitudes has been considerably reduced 
by the fitting. For the EMR the RMS deviation in the 



original models is reduced from 0.56% to 0.20%. For the 
CMR the RMS deviation in the original models is reduced 
even more, from 0.82% to 0.26%. 

Figure [3] shows the evolution of the multipole ratios 
at low Q^ . There is reasonable consistency of the results 
from the different laboratories. The plotted lattice QCD 
results, with a linear pion mass extrapolation[5D] . are in 
general agreement with the data for the EMR but dis- 
agree for the CMR by a wide margin. The EFT analy- 
sis of Pascalutsa and Vanderhaeghen(PV) [H] indicates 
that a linear extrapolation is close to the data for the 
EMR but results in a considerable underestimation of the 
CMR. The results of the two chiral calculations [20ll2T] 
are also presented in Fig. [Sj As these are effective field 
theories they contain empirical low energy constants. For 
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Gail and Hemmert this includes fits to the dominant M^'^ 
multipole for Q^ < o.2 (GeV/c)^ and for the EMR at 
the photon point (Q^ = 0). In order to achieve the good 
overall agreement they had to employ one higher order 
term with another empirical constant. As also can be seen 
from the large estimated errors of the Pascalutsa and Van- 
derhaeghen EFT calculation [21] a treatment of the next 
higher order term is required. In Fig. [3] two representative 
constituent quark models, the newer hypercentral quark 
model (IIQM)|3ij! and an older non-relativistic calcula- 
tion of Capstick [TO], have been included (the relativis- 
tic calculations are in even worse agreement with exper- 
iment). These curves are representative of quark models 

3/2 

which typically under-predict the dominant M-^{ multi- 
pole by ~ 30% and underestimate the EMR and CMR by 
an order of magnitude, even predicting the wrong sign. 
One solution to this problem has been to add pionic de- 
grees of freedom to quark models [17] . All of these models 
treat the Z\ as a bound state and therefore do not have the 
ttN continuum (i.e., no background amplitudes) so that 
cross sections are not calculated. The Sato-Lee [18] and 
DMT [19] dynamical reaction models with pion cloud ef- 
fects bridge this gap and are in qualitative agreement with 
the Q"^ evolution of the data. These models calculate the 
virtual pion cloud contribution dynamically but have an 
empirical parameterization of the inner (quark) core con- 
tribution which gives them some flexibility in these ob- 
servables. By contrast the empirical MAID [23 and SAID 
[29] represent fits to other data with a smooth Q^ depen- 
dence. 

One way to see the major role played by the pion cloud 
contribution to the resonant multipoles is that for this case 
the expected scale for the Q^ evolution is rn^ = 0.02 GeV^. 
In these units the range of the present experiment from Q- 
from 0.060 to 0.20 GeVVc^ is 3 to 10 units. Therefore it is 
not surprising that one should see relatively large changes 
in the predicted Q" evolution of the resonant multipoles as 
is shown in Figure[31 It is also clear that there is significant 
model dependence in these predictions. 

In conclusion, the new data are at the lowest measured 
Q^ for modern clectroproduction where the dominant pi- 
onic contribution is predicted to be increasing. This Q^ 
region is sufficiently low to be able to test chiral effec- 
tive calculations. These results are in qualitative agree- 
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Fig. 1. (Color online) The measured ao = ar + tOh and olt differential cross sections as a function of S^^ at W = 1221 MeV 
and Q^ — 0.060 (GeV/c)^ before (top panels) and after (bottom panels) fitting. The T symbols are our data points and include 
the experimental and model errors (see Table [1} added in quadrature. The EFT predictions [21] are plotted with their estimated 
uncertainties. The other curves represent predictions from the MAID 2003 [28], SL(Sato-Lee)[ig, DMT [19], and SAID [29] 
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models. The lines with dots are the fitted models with the Ej_j_ and Sij_ quadrupole terms set to zero. 
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Fig. 2. (Color online) The measured gtt and ct^t' differential cross sections as a function of 9^^ at W = 1221 MeV and 
Q^ — 0.060 (GeV/c)^. The T symbols are our data points and include the experimental and model errors (see Table [TJ added 
in quadrature. The EFT predictions [21] are plotted with their estimated uncertainties. The other curves represent predictions 
from the MAID 2003 [21], SL(Sato-Lee)[Il], DMT [H], and SAID [21] models. The model curves after fitting are almost identical 
to those before and so have been suppressed. See text for details. 
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merit with lattice calculations with a chiral extrapolation 29. R.A. Arndt et al., Phys. Rev. C66, 055213 (2002), 

to the physical pion mass ^T] , with recent chiral perturba- http : //gwdac . phys . gwu . edu. 

tion theory calculations pOlET] and with dynamical mod- 30. C. Alexandrou et al., Phys. Rev. Lett. 94, 021601 (2005). 

els which explicitly include the pion cloud [l8l|T9]. How- 31. M. De Sanctis et al., Nucl. Phys. A755, 294 (2005). 

ever, all of these calculations require refinements in order 

to obtain quantitative agreement with experiment. This 

includes lattice calculations with lighter pion masses and 

the next order in effective field theory. 
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Table 1. Values of EMR(%), CMR(%), and A/j^ (in IQ^''' /m^+) extracted from these data with three resonant parameter 
fits using the SAID g^, MAID 2003 [11], Sato-Lee(SL) [ig, and DMT [19] models at the A resonance, W = 1232 MeV, at 
Q^ = 0.060 (GeV/c)^. The original model predictions are in the square brackets. The first error is the statistical and cut errors 
added in quadrature. For Af^^ the second error is the systematic error (for the EMR and CMR they are negligible). For the 
average, the third error is the model error defined as the RMS deviation of the results from the four different models. The 
bottom two lines are the EFT predictions of Gail and Hemmert (GH) [20] and Pascalutsa and Vanderhaeghen (PV) [21j . 



EMR (%) 



CMR (%) 



M- 



3/2 



SAID -2.18±0.31[-1.80] 

SL -2.26 ± 0.30[-2.98] 

DMT -2.11±0.28[-2.84] 

MAID -2.56±0.27[-2.16] 



-4.87±0.29[-5.30] 
-4.46 ± 0.25[-3.48] 
-4.85 ± 0.26[-5.74] 
-5.07±0.26[-6.51] 



40.81 ± 0.29 ±0.57[40. 72] 
40.20 ± 0.27 ±0.56[41. 28] 
40.78 ± 0.27 ±0.57[40.81] 
39.51 ± 0.26 ±0.57[40.53] 



Avg. 



-2.28 ± 0.29stat + sys ± 0.20n 



-4.81 ± 0.27stat-Hsys ± 0.26model 40.33 ± 0.27stat ± 0.57sys ± 0.61model 



GH 
PV 



-2.66 
-2.88 ±0.70 



-6.06 
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39.75 ± 3.87 
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Fig. 3. (Color onhne) The low Q^ dependence of the EMR and CMR aX W = 1232 MeV for the 7*p ^ A reaction. The 
T symbols are our data points and include the experimental and model errors (see Table [T} added in quadrature. The other 
data are: the photon point Q [1] and ® [5], CLAS D [6], Bates A [11], Eisner 0[14], and Pospischil ffl [12]. The lattice QCD 
calculations with linear pion mass extrapolations are shown as x [30] and also the recent chiral perturbation calculations of 
Pascalutsa and Vanderhaeghen (PV) (see EFT in Fig. [T| [S] and Gail and Hemmert (GH) (black solid line) [20]. The other 
curves represent the same models as in Fig. [l] The HQM (long-dashed line) [31] and Capstick (short-dashed line) [16] quark 
models have been included. 



